Slow dynamics and stress relaxation in a liquid as an elastic medium 
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We propose a new framework to discuss the transition from exponential relaxation in a liquid to the 
regime of slow dynamics. For the purposes of stress relaxation, we show that a liquid can be treated 
as an elastic medium. We discuss that, on lowering the temperature, the feed-forward interaction 
mechanism between local relaxation events becomes operative, and results in slow relaxation. 



At high temperature, a liquid under external pertur- 
bation relaxes exponentially fast: a relaxing quantity 
q(t) oc exp(— t/r), where r is associated with the tran- 
sition over a single activation barrier. This is known as 
Debye relaxation. On lowering the temperature, relax- 
ation becomes qualitatively different. In many cases, it 
is well approximated by the stretched-exponential func- 
tion, q(t) oc exp(— (t/r)"), where < (3 < 1; smaller 
(3 corresponds to slower relaxation. Slow relaxation is 
universally seen in many supercooled liquids, glasses and 
other disordered systems [H 0, 0|. It describes a very 
sluggish dynamics: in the wide data range, it decays as 
a logarithm of time. The transition to slow relaxation 
signals the onset of glass transformation range [l[ , which 
ends with the dynamical arrest at T g , glass transition 
temperature. T g is conventionally defined from the con- 
dition that relaxation time r exceeds the experimental 
time scale of 100-1000 seconds. Slow relaxation can be 
seen below T g as well, although stronger perturbation 
(e.g. large pressure 0]) is needed to induce relaxation. 

The universal character of slow relaxation strongly sug- 
gests that there should be some fundamental slowing 
down mechanism, that kicks in on lowering the temper- 
ature. Since stretched-exponential relaxation (SER) was 
introduced by Kohhausch in 1854 [5j, it has proved dif- 
ficult to rationalize it from the first principles, without 
invoking postulates or assumptions. This has contributed 
to the lack of consensus about its physical origin [(| 0] • 
For a more detailed discussion, the reader is referred to 
papers that review slow relaxation [D, 0, Q . 

Our thinking of the liquids is shaped by the idea that 
unlike elastic solids, they do not support static stresses. 
At the same time, it has been long known that liquids 
are no different to solids in supporting stresses at very 
high frequencies [l|, at times smaller than system's re- 
laxation time. Whether the general idea of liquids being 
able to relax stress can be related to slow relaxation has 
remained unknown. 

In this paper, we propose a new framework to discuss 
the origin of non-exponential dynamics in a liquid. We 
show that, on temperature decrease, the feed-forward in- 
teraction mechanism between local relaxation events be- 
comes operative, resulting in slow relaxation. This result 
is discussed solely on the basis of elastic response of a 
liquid, without invoking any additional assumptions or 
postulates. 




FIG. 1: Orowan's illustration of a concordant local rearrange- 
ment. Solid and dashed lines around the shaded atoms cor- 
respond to initial and final positions of a rearrangement, re- 
spectively. Arrows show the direction of external stress. 



It is interesting to recall the previous discussion about 
the elastic aspects of liquid behaviour. Some time ago, 
Orowan introduced "condordant" local rearrangement 
events [8j. A concordant local rearrangement is accom- 
panied by a strain agreeing in direction with the applied 
external stress, and reduces the energy and local stress 
(see Figure 1). A discordant rearrangement, on the other 
hand, increases the energy and local stress. This has led 
to a general result that stress relaxation by earlier con- 
cordant events leads to the increase of stress on later 
relaxing regions in a system. Orowan used this result for 
a heterogeneous system, in which stresses are relaxed by 
local rearrangements in certain regions of the system only 
Q . Goldstein applied the same argument to a viscous liq- 
uid 0: consider a system under external stress which is 
counterbalanced by stresses supported by local regions. 
When a local rearrangement to a potential minimum, bi- 
ased by the external stress, occurs (a concordant event), 
this local region supports less stress after the event than 
before; therefore, other local regions in the system should 
support more stress after that event than before [§| . Not- 
ing that "... a molecular theory of this process will be 
extremely complicated", Goldstein proposed that "the 
least any model of the flow process must acknowledge is 
that the extra stress must be supported elsewhere" [§]. 
Note that this proposal implies that a supercooled liq- 
uid should be considered as an elastic medium that can 
support stress relaxation and redistribution. 
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We begin the discussion with introducing local relax- 
ation events (LREs) , which a system uses to relax stress 
locally. These activated events have been given different 
names in the literature (e.g. molecular rearrangements, 
"flow events" etc [![). Because each event carries a mi- 
croscopic change of a macroscopic quantity (e.g. volume 
etc), the number of events governs the dynamics of re- 
laxing observables. To illustrate the argument of Orowan 
and Goldstein above, lets consider the current number of 
LREs n(t), induced in a liquid by an external perturba- 
tion. As discussed by Orowan and Goldstein, because 
an external perturbation introduces bias towards concor- 
dant relaxation events, which support less stress after re- 
laxation, later LREs should support more stress in order 
to counterbalance. The increase of stress on a current re- 
laxing region, Ap, is therefore a monotonously increasing 
function of n. 

It has been shown that activation barrier for a LRE V 
is essentially defined by the elastic energy. In the Shov- 
ing model, for example, V is given by the work of the 
clastic force needed to shove aside the surrounding liquid 
in order for a relaxation event to take place [l(J. When 
the additional stress Ap is loaded on a current relaxing 
region, its relaxation becomes more difficult since more 
work is needed to shove aside the surrounding liquid (or 
open the cage in Figure 1) in order to allow for relaxation 
event to take place. Hence the increase of V is given by 
the increase of work needed to overcome the additional 
barrier created by the elastic force due to the additional 
stress Ap: AV = J Apdq, where the integral is taken 
over the reaction path of a LRE. If q a is the characteris- 
tic activation volume [lfj, AV = Apq a , giving 

V(n) = V Q +q a Ap (1) 

where Vq is the initial value of the barrier. 

Because Ap is a monotonously increasing function of 
n, we find that V(n) is also a monotonously increasing 
function of n. This describes the feed-forward interaction 
mechanism between LREs, in that activation barriers in- 
crease for later events. 

It is important to discuss the condition under which the 
feed-forward interaction mechanism is operative. Let t s 
be the time needed for elastic interaction to propagate be- 
tween LREs, d the distance between neighbouring LREs 
of about 10 A, c the speed of sound, and To the oscilla- 
tion period, or inverse of Debye frequency (to = 0.1 ps). 
Because c = a/ To, where a is the interatomic distance, 
we find that t s — Tod/a, and is temperature- independent. 
On the other hand, relaxation time r, which is also the 
time between two consecutive relaxation events [l(, in- 
creases on lowering the temperature. It is easy to see 
that at high enough temperature, when r = To, t s > t is 
always true because dja > 1. In this case, local events 
relax independently of each other, because at high tem- 
perature, the time between the events is shorter than the 



time needed for elastic interaction to propagate between 
them. Because LREs are independent, we obtain the ex- 
pected high-temperature result that relaxation is expo- 
nential. On cooling the system down, a certain temper- 
ature always gives the opposite condition: t s < r. When 
the time between local relaxation events exceeds the time 
of propagation of elastic interaction between the events, 
local relaxation events do not relax independently, but 
"feel" the presence of each other. 

The maximal distance between those LREs which are 
involved in the elastic feed-forward interaction mecha- 
nism is obtained by setting t s — r: 

d = a— (2) 

TO 

where d is now the range of elastic interactions in a liquid. 

At high temperature, when r = to, d is on the order of 
interatomic distances, but quickly increases on lowering 
the temperature. Since t = tq exp(V/kT), we obtain 

d = aexp(J^j (3) 

where V is the activation barrier. For super- Arrhneius 
increase of V, it is easy to show [ll| that d exceeds the 
system size well above T g . 

It follows from the above discussion that when d ex- 
ceeds the distance between neighbouring LRE in a liquid 
on the order of 10 A, the crossover from exponential to 
non-exponential relaxation takes place due to the feed- 
forward interaction mechanism. The crossover tempera- 
ture T c can be calculated from Eq. (3). 

To calculate how V depends on the current number 
of LREs ri, we introduce the dynamic variable n{t), the 
current number of relaxing events in a sphere of radius 
r = d/2. n(t) starts from zero and increases to its final 
value n T , n(t) — * n T as t — ► oo. Lets consider the cur- 
rent LRE about to relax to be in the centre of the sphere 
(see Figure 2). All previous concordant LREs that are 
located within distance r from the centre, participate in 
the feed-forward interaction, increasing stress on the cen- 
tral region and hence the activation barrier for the central 
LRE. Let Api(0) be the reduction of local stress due to 
a remote concordant LRE i. Generally, Api decays with 
distance, hence we denote Api (r) as its value at distance 
r from the centre. In what follows, we assume, for sim- 
plicity, that Apj(O) are constant, Apj(0) = Ap . The 
increase of stress on the central rearranging region, Ap, 
can be calculated by integrating Api(r): 

fd/2 

Ap = p 4irr 2 Api(r)dr (4) 

J do/2 

where do is on the order of the size of a relaxing region 
(central region in Figure 2), do > 3a (see Figure 1) and 
p is the density of LREs, p = 6n/nd 3 . 
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FIG. 2: Illustration of the feed-forward interaction mecha- 
nism between local relaxation events. This mechanism oper- 
ates within the sphere of radius r = d/2. Shaded and open 
circles represent local relaxing regions inside and outside, re- 
spectively, of the interaction sphere. 

It is important to note that Eq. (4) describes the re- 
distribution of local stresses in an elastic medium. Unlike 
a solid, a liquid does not support static stresses. How- 
ever, it supports stresses at high frequencies in a solid- 
like manner These frequencies correspond to times 
smaller than relaxation time r. We introduce the dis- 
tance over which interactions in a liquid are elastic, liquid 
elasticity length d e i. If c is the speed of sound, d e i = ct. 
Because c = a/ro, we find d e i — ar/ro- As follows from 
Eq. (2), d e i — d. We therefore find that the distance in 
which the feed-forward interaction mechanism between 
LREs operates is also the liquid elasticity length. Hence 
we are justified to use d/2 as the upper limit of integra- 
tion in Eq. (4). 

In an elastic medium, stresses decay as Ap(r) cx 1/r 3 
0. Since Ap(r) = Ap at r = do/2, Ap(r) = 
Apo(do/2r ) 3 . A straightforward integration of Eq. (4), 
combined with Eq. (1), gives 

V^Vo + Vi— (5) 

n r 

where V\ = t: /2p r q a Ap^df l \u{d/ do) and p r = 6n r /(ird 3 ) 
is the density of the final number of events in the sphere. 

Eq. (5) describes the feed-forward interaction mecha- 
nism in a liquid at T < T c . 

We are now set to derive slow relaxation. As discussed 
above, at high temperature above T c , the feed-forward 
interaction mechanism is absent, and LREs relax inde- 
pendently of each other. The rate of LREs, dn/dt, is 
proportional to the number of unrelaxed events, (n r — n), 
and the event probability, p = exp(Vo/fcT), where Vo 
is the high-temperature activation barrier. Introducing 
q = n/n T , we write: 

| (-£)(!-.) W 
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FIG. 3: Solid lines are solutions of Eq. (7), dashed lines 
are fits to SER. From left to right, a=0.3, 1.2 and 2.2. SER 
parameters (J), r) are (0.94, 1.12), (0.78, 1.62) and (0.66, 2.51) 
in the same order. 

where t is dimensionless time t — > t/ro. The solu- 
tion of Eq. (6) is the expected high-temperature expo- 
nential relaxation, q = n/n T = 1 — cxp(— t/r), where 
t = T exp(V /kT). 

As temperature is lowered below T c , the feed- forward 
interaction mechanism between consecutive LREs be- 
comes operative as discussed above. Because this mech- 
anism is described by Eq. (5), the event probability p 
becomes dependent on n (or q), and the rate equation is 
(compare with Eq. (6): 

^ = CX P (~^) ( X -9)exp(-aq) (7) 

where a = V\/kT. 

Unlike Eq. (6), Eq. (7) describes slow relaxation due 
to the factor exp(— aq). Non-exponential relaxation be- 
comes exponential at high temperature when a is small. 
Since SER is one of the simplest functions that describes 
slow relaxation, we attempt to fit the solution of Eq. (7) 
to q = 1 — cxp (— (t/r)P). Because the degree of non- 
exponentiality of the solution of Eq. (7) is controlled by 
a only, we set exp(— Vo/kT) = 1. This affects r but not 
(3. In this form, Eq. (7) is a one-parameter description 
of non-exponential relaxation. We solve Eq. (7) numer- 
ically for different values of a, and fit the solution to 
SER using the least-squares method. A reasonably good 
quality of the fits is illustrated in Figure 3 over about 5 
decades of reduced time. Note that smaller temperature 
(larger a) gives smaller (3 (see the legend in Figure 3), 
consistent with experimental results [2, S] • 

It is interesting to note that on lowering the tem- 
perature, we expect to find the crossover of (5 to the 
lower slope. Using Eq. (3), we find V\ oc ln(d/drj) = 
(ln(a/do) + V/kT). This relation remains true as long as, 
on lowering the temperature, d < L, where L is the sys- 
tem size. For propylene carbonate, for example, a simple 
calculation shows that d = L takes place at temperature 
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above T g [ll|. When d > L, the upper limit of integral 
in Eq. (4) is L, and V\ oc ln(I//do), and is temperature- 
independent. Hence at d = L, V\ shows a kink and starts 
to saturate to a constant value. Because (3 decreases with 
a = Vi/kT (see Eq. (7)), we find that d = L should mark 
the crossover of (3 to the lower slope. This is indeed what 
is seen in several glass-forming systems at low tempera- 
ture 
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We also note that the proposed picture of slow relax- 
ation gives rise to the property that has been discussed 
quite intensely recently, namely that relaxation in super- 
cooled liquids and glasses is dynamically heterogeneous 
[l4l ]. Dynamic heterogeneity naturally arises in our pic- 
ture, since activation barriers increase for later LREs that 
are located in different regions of the system. This gives 
different waiting times in different regions of the system 
and dynamic heterogeneity. 

In this paper, our main intention has been to identify 
and discuss the origin of slow relaxation in a liquid. We 
proposed that slow relaxation can be understood solely 
on the basis of liquid elasticity, without invoking any ad- 
ditional assumptions or postulates. We have shown that 
on lowering the temperature, the feed-forward interac- 
tion mechanism between LREs becomes operative and 
sets slow relaxation. We have offered a formalism to de- 
scribe this mechanism. The mathematical treatment can 
be developed further in future work. Here we note that 
it is an encouraging sign that we can discuss slow relax- 
ation in a liquid using only concepts from elasticity, one 
of the most robust areas of physics. 

Compared with other models of slow relaxation, our 
picture is probably the closest in its spirit to the Cou- 
pling model (CM) of the glass transition The main 
assumption of the CM model is the postulate that non- 
exponentiality is the result of some sort of cooperativity 
of relaxation that kicks in on lowering the temperature. 
Here, we discussed that the origin of cooperativity is the 



feed- forward interaction mechanism between LREs. Ab- 
sent at high temperature, it becomes operative on low- 
ering the temperature, setting the cooperative character 
of relaxation and slow dynamics. 

We have recently shown that the feed-forward interac- 
tion mechanism, which gives non-exponential relaxation, 
is also responsible for the Vogel-Fulcher-Tamman (VFT) 
law for relaxation time. Here, super- Arrhcnius behaviour 
is associated with increase, on lowering the temperature, 
of the range of the feed- forward interaction ll| . Joining 
this result with the present discussion, we can conclude 
that the two most important open questions in the area 
of glass transition, the origin of SER and the VFT law [l[ , 
can be understood on the basis of elasticity of a liquid. 

I am grateful to V. V. Brazhkin, J. C. Dyre and M. 
Goldstein for discussions, and to EPSRC for support. 
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